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Abstract 

We consider the tunneling of a wave packet through a potential barrier which is cou- 
pled to a nonintegrable classical system and study the interplay of classical chaos and 
dissipation in the tunneling dynamics. We show that chaos- assisted tunneling is further 
enhanced by dissipation, while tunneling is suppressed by dissipation when the classical 
subsystem is regular. 
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I. Introduction 



The systems with mixed quantum-classical description have been the subject of considerable 
interest in recent years [1-9]. The validity of this description essentially rests on whether the 
quantum effects of one subsystem is negligible compared to the others. The well-known exam- 
ple is the Maxwell-Bloch equations [3, 4] which describe a model two-level system (quantum 
mechanical) interacting with a strong single mode classical electromagnetic field. The others 
comprise the models involving nuclear collective motion [5] , one dimensional molecule [6] where 
the motion of an electron is described by an effective potential provided by the nuclei and other 
electrons within adiabatic approximation scheme, etc. The mixed description has also been 
employed by Bonilla and Guinea [7] for the measurement processes. Taking recourse to an 
average description in terms of an effective classical Hamiltonian, Pattanayak and Schieve [8] 
have studied some interesting aspects of semi-quantal chaos. 

Very recently the interplay of classical and quantum degrees of freedom in a special class 
of systems with system-operators pertaining to a closed Lie algebra with respect to system 
Hamiltonian has been demonstrated to have acquired special relevance in connection with 
dissipation in quantum evolution [9]. While the quantum degree of freedom is responsible for 
the generic quantum features the implication of classicality is two-fold. First, it has been shown 
that if the classical degree of freedom is assigned the type of evolution as prescribed by the 
classical treatment of dissipativc systems it is possible to realize dissipation for such composite 
systems via an indirect route without any violation of quantum rule. Second, if the classical 
subsystem is nonintegrable then the quantum motion by virtue of this nonintegrability may be 
profoundly affected. Our object here is to study this dissipative evolution of a quantum system 
coupled to a nonintegrable classical system. 

Before going into the further details let us elaborate this issue a httle further in somewhat 
general terms. 

We first consider the coupling of a quantum system to a classical system in terms of the 
following Hamiltonian 

H — Hq + Hci + Hq^cl , (1) 

where Hq and Hd refer to quantum and classical subsystems of the total Hamiltonian, respec- 
tively. Hq-ci is the interaction potential which contains both classical and quantum degrees of 
freedom. If the quantum operators { -Rj | form a closed algebra with respect to the Hamiltonian 
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H of the system, i. e. , if we have relations of the type 

n 

'H{t),Ri\^ihY,9ji{t)Rj i = l,2...n, (2) 

i=i 

where gji are the elements of a n ® n matrix, then one can have, by virtue of the equations of 
motion 

§ = ![//. A] .= 1,2...., (3) 

a set of linear first order differential equations 

^--tgM). (4) 

which describes the temporal evolution of the expectation values. 

It is important to note that Eq. (4) depends crucially on the elements gji which are deter- 
mined by the classical dynamical variables as appeared in the Hamiltonian H (Eq.l). In other 
words, gji is dependent on classical subsystem Hd- If the energy is taken to be the quantum 
expectation value of the Hamiltonian H, then one can easily generate the temporal evolution 
of classical conjugate variables, coordinate x and momentum p using {H) as follows; 

dx _ d{H) 
dt ~ dp ' 
dp d{H) 



dt dx 



(6) 



Eqs. (4-6) comprise the complete dynamical evolution corresponding to Hamiltonian H in the 
dissipation-free case. 

We now take into consideration two different aspects of classical motion of the subsystem 

First, the classical dynamical variable is made dissipative in an ad hoc fashion by adding 
— 7P to Eq.(5), i. e. , 

dp d{H) 

This allows [9] an indirect route to dissipation in the quantum evolution because of Hg-d term 
in Eq.(l). It is easy to check that this quantum description is quite valid since no quantum rule 
is violated in the process. Thus Eqs. 4, 5 and 7 govern the disspative dynamical evolution. We 
note, in passing, that this method has a distinct advantage over some other approaches [10- 
12] to quantum theory of damping which relies on explicit time-dependent Hamiltonian, such 
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as, Kanai Hamiltonian [10] etc., since these approaches lead to contradiction with uncertainty 
principle. 

The second aspect of classicality of the subsystem concerns the nature of motion generated 
by Hci- If the classical subsystem is non-integrable then the quantum evolution is affected 
through the properties of gji matrix. It is important to emphasize a pertinent point at this 
stage. There are well-known cases where exponential instability occurs in a system where the 
overall dynamical system is non-integrable. For example, we refer to the case of optical chaos 
described by Maxwell-Bloch equations [3, 4]. Similar consideration of quantum-classical mixed 
description has been given to explore quantum chaos within Born-Oppenheimer approximation 
in a typical molecular system where the fast electronic motion and the slow nuclear motion 
(classical) separate in a very natural way. In contrast to these cases the present paper concerns 
the nonintegrability due to the subsystem Hd only. 

The studies of dissipative effects in quantum systems are usually based on a traditional 
system-reservoir linear coupling model. There are two standard ways of formulating the 
problem[13-15, 21, 25]. The first one is based on density operator approach and leads to 
the socalled master equation while the second one is based on the Heisenberg picture and leads 
to the noise operators. In the latter case one replaces the reservoir by damping terms in the 
Heisenberg equations of motion for a dissipation-free system and adds fluctuating forces as the 
driving terms which impart fluctuations to the system. The operator forces are such that, flrst, 
the system has the correct statistical properties to agree in the classical limit and second, they 
maintain the commutation properties of the Boson operators to ensure that uncertainty prin- 
ciple is not violated. These consideration are taken care of in the treatments of density based 
or operator based theories (e. g. , see Leggett and co-workers [13] and Krive and coworkers 
[14] and others [15]. The spiritual root of quantum statistical approach to damping lies in the 
fluctuation-dissipation theorem, which illustrates a dynamical balance of inward flow of energy 
due to fluctuations from the reservoir into the system and the outward flow of energy from the 
system to the reservoir due to dissipation of the system mode. Such a dynamical balance is 
automatically maintained in the present treatment (and one need not add a stochastic term 
in Eq.(7)) through a feedback from quantum to classical subsystem since here one deals with 
a finite number of degrees of freedom in the quantum plus classical subsystems as compared 
to the infinite number of degrees of freedom of the reservoir in the traditional system-reservoir 
description. The closure property of the algebra pertaining to the quantum system thus plays 
an important role in the feedback mechanism. Apart from being simple to implement a decisive 
advantage of the approach is that it may be carried over to the nonlinear systems (e. g. , a 
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Morse oscillator described by SU(2) or SU(1, 1) Lie algebra [16]) in a straight forward manner, 
while the master equations are strictly valid for the linear models. 

The consideration of the two aspects of classicality as mentioned earher therefore leads us 
to dissipative chaotic evolution of classical degrees of freedom. Our object is to explore the 
influence of classical chaos on quantum evolution, dissipation being realized indirectly in the 
overall dynamics through classical friction. The generic quantum feature that we study here is 
tunneling. It is extremely important to asses what influence if any classical chaos have on it, 
particularly in presence of dissipation. Although the interplay of chaos and tunneling [17-19] 
or tunneling and dissipation [14, 20, 21] have been the subject of many researchers separately 
over the last decade, it is not clear how dissipative tunneling is influenced by classical chaos. 
We are therefore concerned here with all three interplaying aspects of evolution, e. g. , classical 
chaos, dissipation and tunneling in a typical model with mixed quantum-classical description. 



II. The model and the mixed quantum- classical dynamics 

We consider a particle with fixed energy Eg < penetrating an inverted potential barrier. —Eg 
is the energy measured from the top of the barrier (Fig. 1). The quantum subsystem Hg which 
describes the penetration of the particle through the barrier is given by 

Xg and Pg are the quantum mechanical operators corresponding to position and momentum of 
the particle respectively, m is the mass of the particle and uq refers to the frequency of the 
inverted well. 

The Hamiltonian for the classical subsystem is given by 

Hd = ^ + axil - b^d + 9^d coscjit , (9) 

which governs the motion of a classical system with mass M and characterized by its position Xd 
and momentum variable Pd in a double- well potential driven by an external field with frequency 
LOi. a ad. b are the parameters of the double- well oscillator, g includes the effect of couphng of 
the external field to the oscillator and the strength of the field. 

Hci is nonintegrable and has been widely employed by various workers [17-19, 22-23] over 
the last few years in a variety of situations related to classical and quantum chaos. 
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The description of the Hamiltonian with mixed degrees of freedom is now made complete by 
considering the couphng of classical and quantum degrees of freedom in terms of the interaction 
potential Hg^cii— ^XciXg) so that we have 



H = Hq + Hcl + XXclXg , 

A represents the coupling constant. 

Making use of the following rescaled dimensionless quantities 



(10) 



X — 



(mMwwo)^ 



{Muj)2xa , B = 



Ps 



{Muj)-2p^ , G 



(Ma)3)5 



we rewrite the Hamiltonian(lO) 



"2 "2 2 

t't P ^P A 4 T~l 2 ^ 

H = uiq— — f^Q-^r H — 7: — H f^^s — Boos + Guus cosuJit + x^^o^^ > 



;i2) 



where uj is the linearized frequency of the double-well potential. 

By choosing the relevant operators belonging to the set {i, a;^,^^, L = xp ^ px] we 
may construct a partial Lie algebra (1 is the unity operator). The temporal evolution of 
the expectation values of these operators (see Eq.4) is given by the following set of coupled 
differential equations 



d{x) 

<m 

dr 
d{x^) 

dr 
d{p^) 

dr 
d{L) 

dr 



(L) - 2xs(S>) , 
2((p2^+ {x^)-xs{x)) 



(13) 



where r denotes the scaled dimensionless time r = cuot. 
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Since (H) coincides with the energy of the system, the classical equations of motion are 



ds 

d^ 
dps 
dr 



X{x) + on cos l^iT + AA^s^ - 2Bns + Tp, 



(14) 



where, 



and 



n = — 



r 



UJo 



(15) 



Here F as defined above is the rescaled dimensionless damping constant introduced in the 
classical equations of motion (14) in an ad hoc fashion. We have already pointed out that this 
ad hoc introduction of classical friction leads to dissipation in the overall dynamics without any 
contradiction to uncertainty principle. Eqs. (13) and (14) thus govern the complete dynamics 
of the mixed system. 

We now turn to the motion of the wave packet. In order that the motion of a wave packet 
comes close to the motion of a classical particle it is necessary that its average position and 
momentum follow the laws of classical mechanics. However, this condition is automatically 
satisfied for the inverted harmonic potential barrier we consider here. 

We now describe the wave function of a particle by a Gaussian wave packet of the form 



ip{x, t) = N{t) exp -P{t) {x - a{t)) 



(16) 



where a(t) and (3(t) arc two complex, time-dependent parameters to be determined. A^(t) is a 
normalization factor which is not of much interest here. The Gaussian wave packets have the 
decisive advantage since the ansatz (16) is a solution of the Schrodinger equation. 



H 



4) ^^^''^ 







(17) 



if a and /3 satisfy the following two equations 



and 



^da 



X 



(18) 
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or their scaled version (using r = coot) 

. ^da s(t) a , , 

It is easy to express the expectation values of and others operators in terms of a and 
as follows; 

a/? + OL*P* 



(x) 



+ ' 
2ipp*{a - a*) 

2(3f3 



(3 + (3* 



(L) = ^ . (20) 

The wave function of Gaussian form (16) satisfies the minimum uncertainty condition 

ApAx = ^ . (21) 

with Ax^ — {x^) — {xf and Ap^ — (p^) — {pf . For such a wave packet it is easy to see that 

= . (22) 

Since for a parabolic barrier the energy of classical particle is equal to the average energy of the 
quantum particle one can show that the energy spread for the coupling-free parabolic barrier is 

= a;o' [(P)' V + (^)'Ax'] • (23) 

It is required that the minimum wave packet should have a minimum energy spread Aif| 
also. Then we infer from Eq.(23) that the wave packet must be prepared at the classical turning 
point — xo, where {p) = and (x)^ has a minimum. Let us therefore assume that the center of 
the wave packet reaches the classical turning point at t = with (p) — pq — 0. We then have 
[from Hq^u^-^, Eq.l2] 

euo = ^ (24) 
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so that xo — ±(2e)2, where e{— ^) is a dimensionless parameter (see Fig 1) denoting the 
energy of the particle measured from the top of the barrier. We obtain 

a{t = 0) = ao = —^0 = ^(2e)^ 

(3{t = 0) = f5o = \ . (25) 

The initial conditions (25) imply that the particle's wave packet is as compact as possible when it 
arrives at the turning point. By choosing /? = |, the wave packet satisfies the quantum-classical 
correspondence that the energy of the quantum particle is equal to that of the classical particle 
for parabolic barrier. 

The above two initial conditions for a and /3 also set the initial conditions for the others. 
Thus we further have from Eqs. (20), 

{x)t=o = -(2e)^ , 
{P)t=o = 0, 
{x\=o = ^(l + 4e), 

(P^)t=o = I > 

{L)t=o = . (26) 

Eqs. (26) suggest that the initial conditions required for the dynamics can be manipulated 
by controlling a single parameter e, which denotes the average energy of the wave packet with 
which the particle impinges on the barrier at the left turning point. 

Since the classical subsystem is coupled to the quantum degrees of freedom it is also nec- 
essary to specify the initial conditions for the classical variables for its evolution. To this end 
we fix the parameter set uj — 6.32, A — 0.002, B — 0.25 and G — 0.63 for the present 
study and choose the initial conditions for s and Ps for the regular and chaotic trajectories 
as s = -4.02,-4.52,-5.02 for three regular and s = -8.80,-9.30,-11.31 for three chaotic 
trajectories, Ps being chosen to be zero always. We refer to Lin and Ballentine [13] for a typical 
illustration of the phase space. 

The initial conditions (25) and (26) along with those for the classical variables then allow 
us to follow the evolution of the expectation values and of the wave packet in terms of Eqs. 
(13-14) and Eq. 19. The relevant physical quantities of interest are determined in the next 
section. 

Before closing this section a relevant point need to be discussed here. We have considered 
the regular and chaotic trajectories of the classical subsystem and referred to its phase space as 
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if this subsystem were not coupled to the quantum subsystem. Since the quantum subsystem 
considered here has a potential which is not bounded from below (a prototype model for barrier 
penetration) , the phase space of the mixed system ( quantum plus classical) is unbounded and 
therefore the very notion of classical chaos is far from clear in such a situation. Based on this 
consideration we have selected the chaotic and regular parts of the phase space corresponding 
to the classical subsystem only. 



III. The tunneling probability and the current 

We are now in a position to calculate the wave function and the corresponding probability 
that the incident particle penetrates beyond the position x — (. We write 

T(x>C,r)= rdx\^P{x,T)f . (27) 
Normalization of the wave function(16) leads us to 

\^P{x,T)f = ^^ eM-<r){x-{x)f] , (28) 

where a(T) ^ (3(t) + (3*(r). 

Since a portion of the wave packet has already tunneled through the barrier at r = 0, the 
tunneling contribution due to it should be subtracted from the tunneling probability given by 
(27). Thus between r — and t — oo the wave packet tunnels with the probability 

p — T{x > xo, r — > oo) — T{x > xo,t = 0) 

= 752)pexp[-a(0),^ld.-V^r-"~exp[-a(ocyi<J. • (29) 

V TT JO V TT JO 

Here (x) at r = is defined by {x)q — Xq and (x) at r = oo by (x)^. 

Another important quantity of interest is the tunneling current. The time evolution of the 
tunneling current 

JM = lir'Ji-&'^] , (30) 

2i ox ox 

can be calculated by making use of the wavcpackct(16) after normalization, together with the 
solution for a{t) and f3{t) in terms of Eqs(18) and (19). It is easy to see that the particle reaches 
the end of the tunnel at x = xq where it produces the following current: 

j{xo,r)^-[f3a-{Pay-xo{f3-P*)Mxo,r)f . (31) 
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We would now like to discuss here the following questions ; 

(i) how the tunneling probability (29) depends on classical chaos due to the nonintegrability 
of the classical subsystem. 

(ii) how long it takes a particle to tunnel through the barrier and how it depends on classical 
chaos. 

(iii) since the present model incorporates dissipation through classical friction, it is of interest 
to consider how the dissipative tunneling probability and the current depend on the nature of 
the classical motion of the subsystem. 

Let us first discuss the case of dissipation- free tunneling (F = 0). Fig. 2 depicts the tunnehng 
probability for two sets of regular and chaotic trajectories, when the incident wave packet carries 
the average energy e — 1 . It is evident that irrespective of the initial energy of the classical 
subsystem, the tunneling probability is substantially higher for the chaotic trajectories than 
that for the regular one. The calculation is repeated for lower incident energy (i. e. higher e) 
of the wave packet (e = 3) and the result is shown in Fig. 3. The difference is more marked in 
the region of higher values of coupling constant. 

The tunneling current (31) is shown in Fig. 3 as a function of time for two trajectories 
(one regular and the other chaotic). For a relative comparison we have also plotted the current 
for the coupling-free case. We observe that a particle needs more time to tunnel through the 
barrier when the classical subsystem is regular. 

We now turn to the case of dissipative tunneling (F 7^ 0). In Fig. 4 we display the ef- 
fect of dissipation on the tunneling probability when the classical subsystem is chaotic. A set 
of three chaotic trajectories corresponding to three different initial conditions (ps = 0.0; s = 
— 11.31, —9.30, —8.80) which refer to the widely different energies of the undriven oscillator were 
studied. Similarly the initial conditions for the set of three regular trajectories corresponding 
to widely seperated turning points of the undriven well {ps — 0.0; s — —5.02, —4.52, —4.02) 
were examined. A variation of F was also carried out. For the sake of brevity we have plotted 
here only the representative variation for one chaotic (s = — 11.31, Ps = 0.0) and one regular 
(s = — 5.02, Ps — 0.0) trajectory with and without dissipation. It is interesting to observe that 
dissipation increases the tunneling probability quite significantly when the classical subsystem 
is chaotic (solid curves) but for the regular classical subsystem (dotted curves) tunnehng prob- 
ability is suppressed by dissipation. Such a differential behavior of dissipative tunneling of 
wave packets for the chaotic and regular trajectories (which remains qualitatively the same for 
other sets) is also apparent in the peak height of the tunneling current (Fig. 5), although the 
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tunneling time does not differ too mudi in tfie two cases. 

Summarizing tfie above numerical results we observe that (i) classical chaos of the subsystem 
increases the tunneling probability and decreases the tunneling time quite significantly, (ii) 
Dissipation enhances the tunneling when the classical subsystem is chaotic in contrast to the 
case when the subsystem behaves regularly. 

The problem of chaos-assisted tunneling has been addressed by various workers over many 
years. For example, while studying the model two-dimensional autonomous systems it has been 
observed that energy splitting can increase dramatically with chaos of the intervening chaotic 
layer in which tunneling takes place between distinct, but symmetry related regular phase space 
regimes separated by a chaotic layer [24]. Lin and Ballentine [17] carried out a numerical study 
on a driven double- well oscillator to show that as the separating phase space layer grows more 
chaotic with increasing driving strength the tunneling rate is enhanced by orders of magnitude 
over the rate of the undriven system. Utermann et. al. [15] also investigated the same system 
to point out the role of classical chaotic diffusion as a mediator for barrier tunneling. 

The enhancement of tunneling of wave packet by classical chaos in the present model of 
mixed quantum-classical description can be understood in the light of the of this classical 
chaotic diffusion. Instead of considering tunneling between two regions(tori) mediated by a 
chaotic transport between them we consider here a single tunneling process which starts at one 
of the turning points of the barrier. As the wave packet evolves through the barrier its mean 
motion by virtue of the couphng of the inverted potential with classical subsystem is affected by 
classical chaotic diffusive motion of the latter. While in the former cases tunneling is reversible, 
the process considered here is irreversible in nature. 

The role of dissipation is somewhat more intriguing since one is concerned here with three 
interplaying aspects of evolution, e. g, tunneling, classical chaos and dissipation. The mecha- 
nism of enhancement of tunneling probability by dissipation when the subsystem is chaotic, can 
be understood qualitatively in the following way; It has been shown recently [19] that tunneling 
may get enhanced due to virtual mixing of excited states by dissipative interaction in contrast 
to ground state tunnehng which is suppressed due to dissipation. The more relevant to the 
present work is an earher observation by Ford, Lewis and O'Connell [25] who had employed 
a quantum Langevin equation and found an increased tunneling rate for a particle tunneling 
through a parabolic barrier in a black-body radiation(incoherent) field which behaves as a stan- 
dard harmonic bath. It is important to realize that in the present model the classical chaotic 
subsystem although a few degree-of-freedom system mimics the behaviour of a typical heat bath 
[22] which in course of energy exchange with the quantum subsystem acquires partial quantum 
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character and brings in decoherence in the dynamics. The classical chaotic subsystem is thus 
reminiscent of the background of a black-body radiation field reservoir and the enhancement 
of the resulting tunneling of the wave packet through a parabolic barrier can be understood 
quahtatively in the spirit of Ford, Lewis and O'Connell [25]. 



IV. Conclusion 



In this study we have considered the tunneling of a Gaussian wave packet through a po- 
tential barrier which in turn is coupled to a nonintegrable classical system. The operators 
pertaining to this system with mixed quantum-classical description close a partial Lie algebra 
with respect to the Hamiltonian operator of the system. By introducing a phcnomcnological 
classical friction one realizes a mechanism of dissipation in the overall dynamical evolution 
in this model without violating any quantum rule. Because of nonintegrability the classical 
subsystem admits of chaotic behavior. Wc have studied the interplay of classical chaos and 
dissipation in the tunneling of wave packet through the barrier and shown that chaos-assisted 
tunneling is further enhanced by dissipation while tunneling is suppressed by dissipation when 
the subsystem behaves regularly. Dissipation thus plays a significant role in the evolution of a 
tunneling process in presence of classical chaos. 
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Figure Captions 



Fig. 1. A gaussian wave packet is shown at the left classical turning point —Xq, of an inverted 
parabolic barrier. 

Fig. 2. The tunnehng probabihty(p) is plotted as a function of couphng strength (x), for 
the wave packet's average energy Eg — hu!o{e — 3.0) for different initial positions of classical 
subsystem. The intial positions for three classical chaotic trajectories are (a) s=-11.31, (6) 
s=-9.30, (c) s=-8.80 and for three regular trajectories and (d) s=-5.02, (e) s=-4.52 and (/) 
s=-4.02; Ps — 0. (The quantities are dimensionless; scale arbitrary) 

Fig. 3. The tunneling current (T) is plotted as a function of time(T) for the wave packet's 
average energy e = 1.0, (a) for classical chaotic subsystem, (6) for classical regular subsystem 
and (c) uncoupled system (x)=0. 

Fig. 4. The tunneling probability(p) is plotted as a function coupling strength(x) for the 
wave packet's average energy e = 3 for different damping values (o) r=2.0 and (b) r=0.0 when 
the classical subsystem is chaotic. Similarly curves (c) r=2.0 and (d) r= 0.0 are plotted for 
the regular classical subsystem. 

Fig. 5. The tunneling currcnt(T) is plotted as a function time (r) for the wave packet's 
average energy e = 3 for different damping values (a) r=2.0 and (6) r=0.0 when the classical 
subsystem is chaotic. Similarly curves (c) r=2.0 and (d) r= 0.0 are plotted for the regular 
classical subsystem. 
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